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摘要 : 基于 区 域 的 几何 活动 轮廓 〈 Chan-Vese，CV ) 模型 是 乳腺 超声 图 像 常 用 的 一 种 分 割 算 法 。 但 传统 的 CV 模型 不 能 满 
足 乳 腺 超声 图 像 分 割 精度 高 、 速 度 快 的 要 求 。 因 此 ， 文 章 提出 了 一 种 基于 指数 加 权 平 均 比 率 (Ratio of Exponential 
Weighted Average， ROEWA ) 算 子 改进 的 ”CV 模型 ， 用 于 乳腺 超声 图 像 中 病灶 区 域 的 分 割 。 首 先 ， 计 算 乳 腺 超声 图 像 的 
ROEWA 算 子 。 其 次 ， 基 于 图 像 的 ”ROEWA 算 子 构建 边缘 指示 函数 ， 用 于 代替 CV 模型 中 的 ”Dirac 项 。 最 后 ， 去 除 平滑 项 ， 
从 而 提高 曲线 演化 的 速度 。 实 验 结果 表明 ， 文 章 提 出 的 算法 不 仅 能 提高 分 割 的 精度 ， 而 且 能 显著 提高 分 割 的 速度 。 
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Segmentation of Breast Ultrasound Image with an Improved CV model 
Based on ROE WA Operator 


LIU Lei WEN Tiexiang QIN Wenjian LILing GU Jia 
(Shenzhen Institutes ofAdvanced Technology, Chinese Academy of Sciences, Shenzhen $518055, China) 


Abstract : The Chan-Vese (CV) model has been widely investigated for breast ultrasound (BUS) image segmentation. However, 
the traditional CV model can not meet the requirement of high precision and Speed for BUS segmentation. To address this issue, an 
improved CV model based on the ratio of exponentially weighted averages (ROEWA) operator was proposed in this paper. Firstly, the 
ROEWA of the BUS image was calculated. Then, a ROEWA-based edge indicator function was built to replace the Dirac term of 
traditional CV model. Finally, the smoothing term was removed to improve the speed of curve evolution. Experimental results Show 


that the advantages of the proposed model in terms of computational efficiency and accuracy. 
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乳腺 肿瘤 是 当今 社会 中 一 种 常见 的 女性 多 发 疾病 ， 其 发 病 率 及 死亡 率 均 高 于 其 他 女性 恶性 肿瘤 ， 位 
导 首 位 ， 严 重 威胁 着 我 国 女 性 的 健康 证。 研究 表明 ， 乳 腺 癌 在 早期 发 现时 ， 及 早 诊断 和 治疗 能 够 有 效 地 提 
高 乳腺 癌 的 治愈 率 ， 大 大 降低 死亡 率 巴 。 因此， 如 何 及 早 地 发 现 乳腺 中 的 病变 区 域 ， 实 现 乳腺 瘤 的 早 发 现 、 
早 治疗 ， 从 而 提高 乳腺 癌 的 治愈 率 ， 已 成 为 当今 医学 界 中 一 个 非常 重要 的 课题 。 


收 稿 日 期 : 2017-06-30 修 回 日 期 : 2017-08-24 

基金 项 目 : 国家 自然 科学 基金 项 目 (61472411〉; 深圳 市 基础 研究 项 目 (JCYJ20160429174611494) 

作者 简介 : 刘磊 ， 硕 士 ， 研 究 助理 ， 研 究 方向 为 医学 图 像 处 理 与 分 析 ; 温 铁 祥 ， 博 士 ， 副 研究 员 ， 研 究 方向 为 影像 融合 信息 处 理 ; 秦 文 健 ， 博 士 生 ， 
研究 方向 为 多 模 医 学 影像 处 理 和 计算 机 视觉 ， 李 凌 ， 博 士 ， 高 级 工程 师 ， 研 究 方向 为 立体 视觉 引导 治疗 ;， 训 嘉 〈 通 信 作 者 ) ， 博 士 ， 博 士 研 究 生 导 
师 ， 研 究 方 向 为 数字 病理 智能 分 析 技 术 研 究 ，E-mail: jia.gu@siat.ac.cn。 


07 


超声 影像 技术 凭借 其 适应 面 广 、 价格 低廉 、 无 辐射 、 快 捷 及 无 创 等 特点 ， 逐 渐 成 为 乳腺 肿瘤 早期 筛 查 的 
重要 工具 让。 然而 ， 受 成 像 设 备 和 成 像 机 理 的 影响 ， 超 声 图 像 中 往往 存在 大 量 的 斑点 噪声 ， 使 得 图 像 的 信 
品 比 和 对 比 度 较 低 、 边 缘 变 得 模糊 ， 给 超声 图 像 的 判读 带 来 困难 ,此 外 ， 超 声 图 像 的 判读 和 诊断 主要 取决 于 
医生 个 人 的 临床 经 验 和 主观 判断 ， 带 有 很 大 的 主观 性 。 不 同 的 医生 对 于 同一 幅 图 像 可 能 会 做 出 不 同 的 诊断 
结果 ， 同 一 图 像 在 不 同时 期 也 可 能 会 得 出 不 同 的 结论 。 因 此 ， 为 避免 主观 因素 的 影响 ， 提 高 诊断 结果 的 准 
确 性 和 客观 性 ， 越 来 越 多 的 研究 将 计算 机 辅助 诊断 系统 (Computer-Aided Diagnosis，CAD ) 用 到 乳腺 肿 
瘤 的 临床 诊断 之 中 。 肿瘤 区 域 分 割 是 乳腺 肿瘤 CAD 系统 的 重要 组 成 部 分 ， 是 进行 特征 提取 和 肿瘤 良 、 恶 性 
判别 的 重要 依据 ， 也 是 目前 研究 的 热点 和 难点 。 
活动 轮廓 模型 Active Contour Model，ACM) 是 超声 图 像 分 割 最 常用 的 一 种 分 割 方法 。 根据 曲线 表达 

方式 的 不 同 ， 可 以 分 为 两 类 : 基于 参数 的 活动 轮廓 模型 和 基于 几何 的 活动 轮廓 模型 外 。 其 中 ， 基 于 参数 的 

活动 轮廓 模型 (又 称 为 “Snake 模型 ”) 是 由 Kass 等 中 于 20 世纪 80 年 代 末 提出 ， 该 模型 自 提出 以 来 就 受 

到 了 研究 学 者 们 的 广泛 关注 。Snake 模型 使 用 参数 化 的 形式 表示 曲线 或 曲面 的 形变 ， 把 “目标 边界 提取 问 
> 题 ” 转 化 为 “能 量 泛 函 的 优化 问题 ”。 这 虽然 可 以 提供 闭合 光滑 的 边缘 轮廓 线 ， 但 无 法 自 适 应 地 控制 曲线 
拓扑 结构 的 变化 ， 不 能 准确 地 分 割 出 目标 轮廓 中 的 止 陷 区 域 。 基 于 几何 的 活动 轮廓 模型 是 建立 在 曲线 演化 
理论 和 水 平 集 方法 的 基础 之 上 ， 其 原理 是 先 将 平面 闭合 曲线 隐 含 地 表达 为 三 维 曲 面 函数 的 水 平 集 ， 再 通 
过 曲面 的 变化 来 隐 含 的 求解 演化 的 曲线 。 与 基于 参数 的 活动 轮廓 模型 相 比 ， 基 于 几何 的 活动 轮廓 模型 具有 
自动 处 理 轮廓 线 拓扑 变化 的 优势 。 

比较 经 典 的 几何 活动 轮廓 模型 主要 有 : 1988 年 ，Osher 和 Sethiang 提 出 的 基于 变 分 水 平 集 方法 的 几何 
活动 轮廓 模型 〈 又 称 为 “水 平 集 模型 ”) ， 水 平 集 模型 采用 曲线 的 隐 式 方程 代替 显 式 的 表示 方式 ， 从 而 
使 曲线 的 演化 更 加 灵活 ; 1995 年 ，Casellesm 提 出 的 测 地 活动 轮 廊 模型 《Geodesic Active Contour) ; 2001 
年 ，Chan 和 Veseg 基 于 区 域 型 水 平 集 思想 和 Mumford-Shah 泛 函 理论 提出 的 基于 区 域 的 水 平 集 活 动 轮廓 模 
GE 型 (Chan-Vese, CV) ; 2002 年 ，Chan 和 Vese 中 将 所 提 的 分 片 光 滑 和 分 段 常 值 的 CV 模型 做 了 进一步 改进 ， 
:三 。 “推广 到 了 多 相 水 平 集 的 情况 。 
© 传统 的 CV 模型 利用 图 像 全 局 信息 推动 曲线 演化 ， 没 有 包含 图 像 边缘 的 任何 细节 信息 。 因此 ，CV 模型 

对 初始 轮廓 的 位 置 和 噪声 不 是 很 敏感 ， 对 边缘 模糊 或 缺失 的 图 像 也 能 取得 不 错 的 分 割 效果 。 然 而 ， 传 统 
CV 中 Dirac 函数 的 定义 域 比 较 狭 窄 ， 不 能 有 效 地 检测 远离 演化 曲线 的 目标 边缘 ， 且 需要 进行 重新 初始 化 ， 
计算 量 较 大 、 演 化 速度 较 慢 。 
针对 CV 模型 需要 重新 初始 化 、 演 化 效率 低 ， 不 能 满足 医学 图 像 分 制 的 要 求 这 一 缺点 ， 研 究 者 们 提出 
了 不 同 的 改进 方法 ， 主 要 可 以 分 为 两 类 : 一 类 是 在 分 割 的 过 程 中 增加 窗口 化 处 理 W， 另 一 类 是 调整 演化 的 
偏 微分 方程 由 。 其 中 ， 在 分 割 的 过 程 中 进行 窗口 化 处 理 ， 可 以 使 得 活动 轮廓 的 演化 速度 随 着 窗口 的 规模 进 
行 调整 ， 减 少 演变 所 需 的 迭代 次 数 和 计算 量 。 调 整 演化 的 偏 微分 方程 部 分 最 常见 的 做 法 是 :构建 边缘 指示 
函数 g(R) 蔡 换 Dirac 函数 ， 以 消除 Dirac 函数 对 检测 远离 活动 轮廓 线 C 的 抑制 ， 去 掉 影 响 演化 速度 的 平滑 
项 ， 以 提高 演化 的 速度 ， 提 高 分 割 效 率 。 其 中 ， 边 缘 指示 函数 的 构建 大 多 数 是 基于 图 像 的 梯度 信息 2。 

与 以 往 改 进 的 CV 模型 方式 不 同 ， 本 文 提出 了 一 种 基于 指数 加 权 平 均 比 率 (Ratio of Exponential 

Weighted Average， ROEWA) 算 子 改进 的 CV 模型 。 此 处 基于 ROEWA 算 子 而 非 高 斯 边缘 梯度 算 子 的 原因 
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在 于 ， 超 声 图 像 中 固有 的 斑点 噪声 和 合成 孔径 雷达 〈Synthetic Aperture Radar，SAR) 图 像 中 的 噪声 都 是 
乘 性 噪声 ， 而 ROEWA 算 子 较 高 斯 边缘 梯度 算 子 更 适用 于 含有 乘 性 噪声 的 图 像 。 主 要 的 改进 之 处 在 于 : 基 
于 图 像 的 ROEWA 算 子 构建 边缘 指示 函数 ， 用 于 代 蔡 CV 模型 中 的 Dirac 项 ， 改 进 了 CV 模型 的 偏 微分 方 


程 ， 使 分 割 算 法 更 好 地 适应 乳腺 超声 图 像 。 
2 传统 CV 模型 及 ROEWA 算 子 


2.1 传统 CV 模型 

CV 模型 ， 又 称 为 基于 区 域 的 活动 轮廓 模型 ， 是 由 Chan 和 Vese 于 2001 年 提出 的 一 种 基于 图 像 全 局 区 
域 信息 的 能 量 模型 四 。 其 核心 思想 是 利用 分 段 常数 去 逼近 图 像 的 各 个 部 分 ， 通 过 变 分 法 引入 水 平 集 建 立方 
程 ， 采 用 差分 方法 进行 数值 计算 ， 最 终 将 边界 轮廓 的 演化 问题 转化 为 能 量 最 小 化 问题 。 对 于 给 定 的 竺 分 割 
图 像 ”7x,y) ， Chan 和 Vese 定 义 其 能 量 函数 为 ; 


< 


El(c,c,C)=H 虹 ength(C) v Arealinside(O)| 
十 碳 ice CR y) ol dxdy 


Gy) cl dxdy 


2 outside(C ) 


其 中 ，y20 、v20 、H>0、4>0 是 各 个 能 量 项 的 权重 系数 ，c: 和 ec 分 别 表示 活动 轮廓 线 C 内 部 区 域 和 外 
部 区 域 的 平均 灰 度 Length(C) 和 4rea[inside(C)] 分 别 表示 闭合 活动 轮廓 线 C 的 长 度 和 内 部 区 域 的 面积 ， 
用 于 平滑 演化 曲线 。 CV 模型 就 是 寻找 最 优 划 分 C， 从 而 使 得 原 图 像 foc y) 被 活动 轮廓 线 C 划分 为 目标 和 
背景 两 个 区 域 。 故 CV 模型 利用 了 图 像 的 全 局 区 域 信息 ， 通 过 优化 能 量 公 式 E(ci, cs,C) 以 求 得 最 终 的 分 


制 曲 线 。 


在 数值 求解 过 程 中 ， 采 用 水 平 集 函数 C 的 零 水 平 集 表示 活动 轮廓 曲线 C (CeQ， 人 QQ 表示 整 个 图 像 
域 ) ， 则 : 


[x| 


类 一 fx,y) Q19(x,y) 0 
inside(C) 一 fx,y) 师 |9(x,y) 0} 
outside(C)= {x, y) 师 |9(x,y) 0} 
活动 轮廓 线 C 。 内 、 外 区 域 的 平均 灰 度 cl 和 co 的 表达 式 如 下 : 
_ Tox,y)H (Gdxdy 
H(G)dxdy 


1 


_ LC, WUT Hb) dxdy 
[TH(G)]dxdy 


其 中 ， 有 (9) 通常 采用 Heaviside 函 数 ， 其 表达 式 如 下 : 


2 


闻 ( 六 = 二 J et 
2 区 CE 


ht 


其 中 ， 为 常数 。 
故 ， CV 模型 的 能 量 公式 〈 1 ) 可 以 表示 为 : 
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E(c,c,@)=4 ,519Gx, | $Cx,y)dxrdy 
+v H[GG, wldxdy 
+ 7x,y) ol HIGG, y)]axdy 
十 五 直人 ol {1 HlGG, yrdy 


这 样 ， 求 解 上 述 能 量 函数 的 最 小 化 问题 就 转变 成 了 求解 妃 关 于 ce 、c 和 & 的 极 小 值 问 题 。 若 要 求 得 使 能 量 
函数 五 取 极 小 值 时 的 ce 、c: ， 只 需 分 别 求 出 五 关于 cl 、@c 的 一 阶 导数 ， 并 令 其 为 0， 解 出 ce 、c 的 值 即 可 。 
为 了 求解 使 E 取 得 极 小 值 时 的 水 平 集 函 数 包 令 cl 、cz 的 值 不 变 ， 利 用 变 分 法 得 到 其 欧 拉 - 拉 格 衣 日 方程 
引入 时 间 变 量 t 并 采用 标准 的 梯度 下 降 法 进行 迭代 求解 ， 最 终 得 到 控制 水 平 集 函 数 演化 的 偏 微 分 方程 : 


关 -6t) 合 芥 Vv 4 (wu a) 1 (uo c,) 
tk bo (x, 了) 


eal 


其 中 ，65,.(9) 是 五 (9) 的 一 阶 导 数 ， 其 表达 式 为 : 
qd 0, I |e 
SO 
TE ++2 


2.2 ROEWA 算 子 

ROEWA 算 子 是 一 种 基于 线性 最 小 均 方 误差 (Minimum Mean Square Error，MMSE) 的 指数 加 权 均 值 
比率 滤波 器 ， 多 用 于 合成 孔径 雷达 (SAR) 图 像 的 边缘 检测 鸣 ， 由 Fjortoft 等 于 1998 年 提出 。 基 于 线 
性 最 小 均 方 误差 的 指数 滤波 器 在 估计 检测 窗口 内 的 局 部 均值 和 时， 计算 的 是 加 权 后 的 均值 ， 而 不 是 算术 均值 ， 
并 且 这 种 均值 的 变化 与 原 图 像 的 标准 差 成 正比 ， 因 此 ROEWA 检测 具有 恒定 虚 警 率 ( Constant False Alarm 
Rate， CFAR ) 特性 。 
在 一 维 情况 下 ， 该 滤波 器 的 表达 形式 为 : 


fp)=C eH 
其 中 ，C 为 归 一 化 常数 ， w 为 滤波 系数 ，a 的 值 与 噪声 抑制 成 反比 ， 与 边缘 保持 成 正比 。 
在 离散 情况 下 ， 滤 波 器 f(x) 可 以 用 因果 滤波 器 f(x) 和 非 因 果 户 (x) 来 实现 ， 其 离散 形式 下 的 表达 


式 为 : 
b 
f (X= Ef) | TT) 
其 中 ， 


/i(X)=a bh(x) 


六 oO 一 oO hx) 


其 中 ，0<< 括 所 ”1 ; h(x) 是 离散 单位 阶 路 函数 ， 通 常 选用 Heaviside 函 数 ， 当 x>0 时 ，h(x) 三 1， 否 则 


为 0。 
根据 f(x,y) 二 A(X)f(y) 将 上 述 滤 波 器 从 一 维 空间 推广 到 二 维 空间 。 例 如 ， 若 要 计算 图 像 水 平方 向 的 
边缘 强度 分 量 ， 先 使 用 一 维 平滑 滤波 器 了 对 图 像 7(x, y) 逐 列 进行 平滑 滤波 ， 然 后 分 别 用 因果 滤波 器 


-~~ 
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f(x) 和 非 因果 滤波 器 有 (x) 对 平滑 滤波 的 结果 逐 行 进行 滤波 ， 得 到 图 像 水 平方 向 的 无 偏 估 计 Lx (x,») 
和 Lx, (Xx, y) 
Rx (x, yp)=f (xX) *[f (yg (x,y)] 


hy, (x, yp)=f (x) * [f(y)9 (x,y)] 
其 中 ，* 为 水 平方 向 上 的 卷 积 ; 9 为 垂直 方向 上 的 卷 积 。 同 理 ， 将 滤波 顺序 颠倒 可 得 垂直 方向 上 的 无 偏 估 


计 必 (x, yp) 和 boy (Xx, ») 
bs (x, p=f(V) fx) T(x, y)] 


bys (Xs)=f (yf (OO* Tx, y)] 


得 到 两 方向 的 无 偏 估计 后 ， 根 据 ROEWA 算 子 定 义 比值 可 以 求 得 两 方向 上 的 梯度 估计 。 ROEWA 算 子 定义 
比值 为 : 


纹 ? L 
太一 max 开 一 ,一 
让 


于 光学 图 像 的 基于 梯度 的 边缘 检测 算法 ， 可 求 得 2 维 图 像 中 的 总 体 边缘 强度 为 : 


R(x, 人 人) 一眼 CI) Ra (x y) 


3. 基于 ROEWA 算 子 的 改进 CV 模型 


淹 
也 
浴 


针对 传统 CV 模型 中 Dirac 函数 定义 域 较 狭 窄 ， 不 能 有 效 地 检测 远离 演化 曲线 的 目标 边缘 ， 且 需要 重 
新 初始 化 、 计 算 量 较 大 这 一 缺点 ， 本 文 拟 对 其 进行 改进 一 一 构建 基于 图 像 边缘 强度 信息 的 边缘 指示 函数 
&8(R) 取代 传统 CV 模型 中 规则 化 的 Dirac 函数 。 此 处 边缘 函数 旨 在 检测 远离 演化 曲线 的 目标 边缘 ， 引 时 
演化 曲线 停止 在 目标 轮廓 的 边界 处 ， 提 高 收敛 速度 。 基 于 这 个 目的 ， 本 文 构建 的 边缘 指示 函数 如 下 : 


g(R) 一 ec 


其 中 ，wx 是 图 像 边 界 强度 场 对 演化 速度 的 控制 系数 ， 可 取 值 任意 实数 。 在 具体 实验 时 ， 取 值 应 适中 ， 过 大 
的 取 值 容易 使 轮廓 还 未 扩展 到 实际 边界 处 就 已 经 停止 ， 过 小 则 会 使 轮廓 线 穿越 物体 弱 边 界 处 。R 为 由 
ROEWA 算 子 得 到 的 边缘 强度 。 这 里 采用 ROEWA 算 子 而 非 一 般 图 像 分 割 算法 采用 的 梯度 值 ， 主 要 是 考虑 
到 超声 图 像 中 的 斑点 噪声 为 乘 性 噪声 ， 而 ROEWA 算 子 能 更 好 地 适应 噪声 为 乘 性 模型 的 图 像 。 所 构建 的 边 
缘 指示 函数 8(A) 是 一 个 单调 递减 函数 ， 当 R 一 w% 时， g(R) 一 0 ， 演 化 停止 ， 表 明 该 处 越 接 近 于 肿瘤 的 真 
实 边界 轮廓 。 
本 文中 8g(R) 采用 的 是 负 指 数 函 数 形式 ， 而 非 倒 数 形式 ， 是 因为 负 指数 函数 形式 相 比 倒数 形式 具有 更 
好 的 副 零 效果 。 改 进 后 的 CV 模型 的 表达 式 为 : 

剖 一 (a 有 Y A(uo a) (uo c,)] 


$0, x， y) 二 fx, y) 
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SE 


颈 
下 请 项 div 久 - 计算 量 较 大 ， 且 有 对 水 平 集 本 数 有 重新 初 化 的 要 求 ， 严 重 限制 了 本 数 的 演化 加 度 ， 


敛 速度 ， 故 本 文 在 实际 运算 中 忽略 了 该 平滑 项 。 最 终 采 用 的 偏 微分 方程 为 : 


=gtR tht oj Xl a y] 
办 0, x, y)= B(x, y) 


(21) 


4. 实验 结果 与 分 析 


通过 实验 对 比 发 现 ， 去 掉 该 项 后 ， 对 分 割 精度 并 无 影响 ， 因 为 边缘 指示 函数 本 身 也 可 以 加 快 CV 模型 的 收 


实验 所 采用 的 仿真 软件 平台 是 MATLAB R2013a ， 系 统 平台 为 64 位 的 Windows10 操 作 系 统 ，CPU 时 钟 
频率 为 2.94GHz ， 内 存 为 4.0GB 。 实 验 数 据 均 由 合作 医院 广东 广州 中 山大 学 附属 第 三 医院 超声 科 提 供 ， 共 
含 61 例 ， 其 中 恶性 肿瘤 32 例 ， 良 性 肿瘤 29 例 ， 所 有 的 病例 均 有 相应 的 病理 诊断 报告 。 如 无 特殊 说 明 ， 参 数 
设置 如 下 : 时 间 步 长 rz= 0.1 ， 空 间 步 长 1=1、1=4=1、A= 02 、v=0、vwx=2， 友 代 次 数 


n=200 。 


为 了 获得 比较 理想 的 分 割 效 果 ， 所 有 的 图 像 在 分 割 前 均 采 用 各 项 异性 扩散 〈Speckle Reducing 


Anisotropic Diffusion ， SRAD ) 模型 中 进行 了 预 处 理 。 在 本 文 试验 中 ， SRAD 模型 中 各 参数 设置 为 : 
扩散 方程 为 指数 式 ， 迁 代 次 数 寻 二 30， 传 导 系数 上 =30、 和 = 1 。 


为 了 验证 基于 像素 点 的 比率 值 比 一 般 常 用 图 像 分 割 中 常用 的 梯度 值 能 更 好 地 适用 于 超声 图 像 ， 将 图 像 


分 割 中 常用 的 基于 高 斯 边缘 梯度 算 子 的 边缘 指示 函数 与 本 文 基于 ROEWA 算 子 的 边缘 指数 函数 进行 了 比较 ， 


结果 如 所 示 。 从 中 可 看 出 ， 本 文 所 设计 的 边缘 指示 函数 能 更 好 、 更 清晰 地 检测 出 肿 
瘤 的 边界 区 域 。 这 是 因为 超声 图 像 中 携带 的 斑点 噪声 主要 是 乘 性 噪声 ， 而 ROEWA 算 子 对 乘 性 噪声 更 加 适 


( b ) 基于 高 斯 边缘 梯度 算 子 的 边缘 指示 函数 
边缘 指示 函数 
图 1 两 种 不 同 边缘 指示 函数 对 比 图 
Fig. 1 A comparison between two different edge detectors 


为 验证 本 文 算法 的 分 割 精度 ， 采 用 传统 的 CV 模型 、 基 于 高 斯 边缘 梯度 算 子 改进 后 的 模型 和 本 文 提 


的 算法 对 经 病理 确诊 后 的 乳腺 超声 图 像 进行 分 制 ， 并 与 医生 手动 分 割 结果 作 比 较 ， 分 制 结果 如 图 2 所 示 。 


( a ) 初始 轮廓 示意 图 


| 
| 
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( b ) 传统 CV 模型 的 分 割 结果 


( c ) 基于 高 斯 边缘 梯度 算 子 改进 后 的 模型 分 割 结果 


( qd ) 本 文 方法 效果 图 


〈《e ) 医生 手动 分 割 结果 


图 2 ”乳腺 肿瘤 超声 图 像 分 割 结果 对 比 图 
Fig. 2 Comparison of breast tumor segmentation in clinical ultrasound image 

从 图 2 可 看 出 ， 在 给 定 初始 轮廓 的 前 提 下 ， 传 统 的 CV 模型 和 基于 高 斯 边缘 梯度 改进 后 的 模型 均 能 
检测 出 肿瘤 区 域 的 整体 轮廓 ， 但 在 一 些 细节 上 的 处 理 效果 并 不 太 理 想 。 此 外 ， 人 针对 部 分 图 像 ， 随 着 迭代 次 
数 的 增加 ， 演 化 曲线 会 越过 肿瘤 边界 继续 向 外 演化 ， 导 致 过 分 割 而 生成 错误 边界 ， 如 图 2 〈b ) 和 图 3 
所 示 。 而 如 图 2 〈d) 所 示 ， 采 用 本 文 算 法 ， 不 仅 可 以 更 好 地 分 割 出 肿瘤 区 域 的 边界 轮廓 ， 而 且 随 着 迭 
代 次 数 的 增加 ， 边 界 可 以 稳定 地 固定 在 边缘 处 而 不 会 过 分 割 。 

为 了 定量 评价 本 文 方法 对 超声 图 像 分 割 的 质量 ， 以 医生 手动 分 割 的 结果 《〈 图 2 (e) ) 作为 实验 结果 
参考 的 金 标准 。 由 于 肿瘤 区 域内 部 往往 出 现 一 些 钙化 ， 在 图 中 表现 为 高 回声 区 域 ， 而 在 分 割 的 过 程 中 ， 这 
部 分 区 域 往往 会 被 单独 分 割 出 来 ， 表 现 为 一 些小 空洞 ， 如 图 2(b) (d) 所 示 。 在 临床 上 ， 这 些 “ 高 回声 区 
域 ”属于 病变 区 域 ， 故 在 进行 性 能 比较 时 ， 比 较 的 是 空洞 填充 后 的 二 值 化 图 像 。 采 用 准确 度 、 灵 敏 度 、 特 
异性 、 阳 性 预测 值 ( Positive Predictive Value ) 和 阴性 预测 值 ( Negative Predictive Value ) 共 5 
个 指标 来 量化 评价 分 割 的 性 能 ， 结 果 如 表 1 所 示 。 


表 1 几 种 分 割 算法 的 性 能 比较 


Table 1 The performance comparison of several segmentation algorithm 


， 定量 指标 
活 全 准确 度 灵敏 度 特异 性 阳性 预测 值 阴性 预测 值 
传统 CV 模型 0.856 5+0.101 0.891 14+0.1079 0.8086+0.2038 0.6365+0.1307 0.8916+0.2215 
基于 梯度 改进 后 的 模 0.965 3+0.014 0.811 5+0.057 3 0.8962+0.2227 0.9703+0.0185 0.8689+0.2168 
本 文 算法 0.974 2+0.006 0.878 4+0.020 6 0.9909+0.0050 0.9563+0.0045 0.9827+0.0052 


为 验证 本 文 算法 在 分 割 速度 方面 的 提升 ， 采 用 传统 CV 模型 和 本 文 改进 后 的 CV 模型 对 同一 幅 图 像 进 
行 分 割 ， 并 将 处 理 过 程 中 不 同 迭 代 次 数 的 曲线 演化 情况 和 所 用 时 间 进 行 统计 ， 结 果 如 和 所 示 。 从 表 2 可 以 
看 出 ， 与 传统 的 CV 模型 相 比 ， 本 文 改进 的 CV 模型 显著 地 提高 了 曲线 演化 的 速度 ， 可 以 更 快 地 演化 到 肿 
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瘤 的 真实 边界 。 也 就 是 说 ， 本 文 改 进 后 的 CV 模型 可 以 极 大 地 提 


Sa 


高 图 像 分 割 的 效率 。 


ca 


传统 CV 模 型 荐 3 


过 代 ”60 次 过 代 ”100 次 迁 代 200 次 连 代 300 次 
图 3 ”传统 CV 模型 和 本 文 改进 后 CY 模型 的 演化 速度 对 比 图 
Fig.3 Comparison of the standard CV model with the CV-ROEWA model 

表 2 不 同 迭 代 次 数 下 的 分 割 时 间 对 比 


Table 2 Comparison of segmentation time in different iterations 


迁 代 次 数 时 间 ( s) 
传统 CV 模型 本 文 模型 
60 12.672 3.028 
100 19.352 4.258 
200 27.251 8.037 
300 44.208 11.815 


5， 小结 


肿瘤 区 域 分 割 是 医学 图 像 CAD 系统 的 重要 组 成 部 分 ， 精 确 地 提取 病变 区 域 的 边缘 轮 廊 有 利于 准确 提 
取 其 形态 特征 ， 从 而 提高 肿瘤 良 恶性 判别 的 准确 度 。 针 对 乳腺 超声 图 像 的 分 割 难题 ， 本 文 在 充分 分 析 传 统 
CV 模型 优点 和 缺点 的 基础 上 ， 提 出 了 一 种 基于 ROEWA 算 子 的 改进 CV 模型 。 实 验 结果 表明 ， 本 文 算法 
不 仅 能 提高 分 割 的 准确 度 ， 而 且 能 显赫 提高 分 割 的 速度 。 
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